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ABSTRACT 



We report the results of the spectral analysis of two observations of the Vela 
pulsar with the Chandra X-ray observatory. The spectrum of the pulsar does 
not show statistically significant spectral hues in the observed 0.25-8.0 keV 
band. Similar to middle-aged pulsars with detected thermal emission, the 
spectrum consists of two distinct components. The softer component can 
be modeled as a magnetic hydrogen atmosphere spectrum — for the pulsar 
magnetic field B — 3 x 10^^ G and neutron star mass M — 1.4Mq and radius 

= 13 km, we obtain = 0.68 ± 0.03 MK, Lg^i = (2.6 ± 0.2) x 10^2 erg s-\ 
c? = 210 ± 20 pc (the effective temperature, bolometric luminosity, and radius 
are as measured by a distant observer). The effective temperature is lower than 
that predicted by standard neutron star cooling models. A standard blackbody 
fit gives = 1.49 ± 0.04 MK, L^^^ = (1.5 ± 0.4) x 10^2 d^^Q erg s"^ (cisso 
is the distance in units of 250 pc); the blackbody temperature corresponds 
to a radius, R°° — (2.1 ± 0.2) (i25o km, much smaller than realistic neutron 
star radii. The harder component can be modeled as a power-law spectrum, 
with parameters depending on the model adopted for the soft component 
— 7 = 1.5 ± 0.3, = (1.5 ± 0.4) x 10^^ dj^o erg s"^ and 7 = 2.7 ± 0.4, 
Lj. — (4.2 ± 0.6) X 10^^ (^250 erg s"-*^ for the hydrogen atmosphere and blackbody 
soft component, respectively (7 is the photon index, is the luminosity in the 
0.2-8 keV band). The extrapolation of the power-law component of the former 
fit towards lower energies matches the optical flux at 7 ~ 1.35-1.45. 



Subject headings: stars: neutron — pulsars: individual (PSR B0833-45) — 
X-rays: stars 
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1. Introduction 

X-ray observations of rotation-powered pulsars with the Einstein, EXOSAT and, 
particularly, ROSAT observatories have estabhshed (Ogelman 1995) that at least three 
middle-aged pulsars, PSR B0656-I-14, B1055-52 and Geminga, with characteristic ages of 
10^-10^ yr, show thermal soft X-ray radiation, with (blackbody) temperatures 0.5-1 MK, 
interpreted as emitted from the surfaces of cooling neutron stars (NSs). Investigation of 
this radiation is needed to understand the NS cooling history and to study the properties 
of the NS surface layers. Further observations with the ASCA observatory have confirmed 
that, as it had been suggested from the ROSAT data, the spectra of these pulsars at 
higher energies ( > 1.5-2 keV) are dominated by nonthermal, power-law components, with 
photon indices 7 ~ 1.3-1.6, presumably generated in the NS magnetospheres (e.g., Wang et 
al. 1998). Studying the power- law components is important to understand the mechanisms 
of the multiwavelength (optical through gamma-rays) pulsar radiation. 

In addition to these three pulsars (dubbed "The Three Musketeers" by Becker & 
Triimper 1997), the (presumably) much younger Vela pulsar (r = P/[2P] — 11 kyr, 
P = 89.3 ms, E = 6.9 X 10^^ erg s-\ 5 ~ 3 x 10^^ G — see Taylor, Manchester, & Lyne 
1993 for other properties and references) has often been mentioned as a possible source 
of thermal radiation, based on softness of its spectrum (e.g., Ogelman 1995). Ogelman, 
Finley, & Zimmermann (1993) analyzed the spectrum of the Vela pulsar detected with 
ROSAT. They obtained acceptable fits of the spectrum with a power law (PL) model 
with 7 ~ 3.3 and a blackbody (BB) model of temperature T°° ~ 1.7 MK. The slope of 
the spectrum obtained in the PL fit is much steeper than typical nonthermal spectra of 
young and middle-aged pulsars. The temperature and luminosity inferred from the BB fit 
correspond to an effective radius of emitting region about 1.3 km at a distance d — 250 pc 
(Cha, Sembach, & Banks 1999). A two-component fit with a BB-I-PL model yielded a lower 
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temperature, 1.3 MK, a larger radius, 3 km at 250 pc, and a surprisingly small photon 
index, 7 ~ 0.1 (Ogelman 1995). Page, Shibanov, & Zavlin (1996) fit the same spectrum 
with the NS hydrogen atmosphere models (Pavlov et al. 1995) and obtained a lower effective 
temperature and a larger R/d: = QrT^s ^ 0.8 MK, d ~ 300 pc for R°° = g'^R = 13 km 
(qj. = [1 — 2GM / Rc^Yl"^ is the gravitational redshift parameter). 

However, the true spectrum of the observed soft X-ray pulsar radiation (not to mention 
the presence of two components) has remained elusive because angular resolution of the 
X-ray telescopes has been too low to separate the pulsar radiation from that of the bright 
pulsar-wind nebula (PWN) of about 2' diameter around the pulsar (Harnden, Grant, & 
Seward 1985; Markwardt & Ogelman 1988). The subarcsecond resolution of the Chandra. 
X-ray Observatory (CXO) provides the first opportunity to resolve the pulsar spectrum. 
The high energy resolution of the CXO grating spectrometers allows one to address another 
important problem — if the soft component of the pulsar radiation indeed originates 
from the NS surface layers (an atmosphere), it may show spectral features (lines and 
photoionization edges) which can be used to investigate chemical composition, gravity and 
magnetic field of the surface layers. 

In this Letter wc present first results on the X-ray spectrum of the Vela pulsar obtained 
with CXO. We describe the observations and the spectral analysis in §2 and discuss some 
implications of the observed spectrum in §3. 

2. Observations and Spectral Analysis 

To investigate the low-energy part of the pulsar spectrum with high energy resolution, 
the spectrum dispersed with the Low- Energy Transmission Grating (LETG; see Brinkman 
et al. 2000) was imaged on the Spectroscopic Array of the High Resolution Gamera 
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(HRC-S; see Murray et al. 1997). The 25.6 ks observation (ObsID 127) was taken on 28 
January 2000. For the analysis we used the Level 2 data files produced by the standard 
pipeline processing (v. R4CU5UPD11.1). The dispersed LETG spectrum (dispersion 
1.148 A/mm) was extracted from a strip of a 3'.'6 (0.18 mm) width in the cross-dispersion 
direction, which contains about 90% of source events dispersed at a given wavelength 
(see Chandra Proposers' Observatory Guide, v.3.0; CPOG hereafter). The zero-order and 
dispersed radiation from the PWN contaminates the dispersed pulsar spectrum at short 
wavelengths, A < 5-6 A, and the plate background dominates at A > 60 A. Since at 
A = 50-60 A there are data only from the positive dispersed orders, we chose the range 
A = 6-50 A (E = 0.25-2.0 keV) for spectral analysis. The background was taken from 
boxes between 10" and 30" (0.49-1.46 mm) from the center of the source spectrum in the 
cross-dispersion direction. To search for spectral lines in the dispersed pulsar spectrum, 
wc binned the extracted sourcc-plus-background and background spectra in 0.02 A bins 
(intrinsic resolution of the instrument in the chosen wavelength range is ~ 0.05 A). We 
used groups of 10-20 sequent bins to estimate deviation in number of source counts in each 
bin of a given group from the mean value in the group. The maximum deviation, in both 
the positive and negative dispersed orders, was found to be at a 2.7cr significance level. 
Therefore, we conclude that the LETG spectrum of the Vela pulsar shows no statistically 
significant spectral lines. 

For further analysis, we removed contribution of the higher orders at longer wavelengths 
using a "bootstrap" method which works as follows (see also CPOG). One assumes that 
contamination from the higher orders can be neglected at the shortest wavelengths. The 
source counts in a short wavelength bin are used to deduce the corresponding high-order 
contributions of this bin to longer wavelengths with the aid of relative grating efficiencies of 
the higher orders. These contributions are subtracted, and the process is repeated working 
upwards in wavelength. Using this method, we obtained the first-order pulsar spectrum 
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(91% of the total dispersed data), with a source countrate 281 ± 6 ks ^ (312 ± 6 ks ^ after 
correcting for the extraction efficiency) in the 0.25-2.0 keV range 

To examine various continuum models, we grouped the counts in 162 energy bins with 
at least 30 source counts per bin and fit the models making use of the first-order effective 
area0. We found that a simple PL model requires a large photon index 7 = 4.1-4.3 and 
an implausibly high hydrogen column density nH,2o = n-n/lO"^^ cm~^ = 9.1-9.3; moreover, 
the best-fit model spectrum exceeds the observed one at -E > 1.2 keV. The spectrum fits 
much better with thermal models (blackbody, NS hydrogen atmosphere), but the observed 
spectrum somewhat exceeds the model spectra at the highest energy channels (1.5-2.0 
keV), which indicates the presence of a second component with a harder spectrum. 

To search for the harder component, we used the archived 37.0 ks observation^ 
(ObsID 131; pipeline processing v. R4CU5UPD8.2) taken on 11-12 October 1999 with 
the Spectroscopic Array of the Advanced CCD Imaging Spectrometer (ACIS-S; Garmire 
et al. 2001) in Continuos Clocking (CC) mode which allows timing at the expense of 
one dimension of spatial resolution (along the CCD chip columns). This observation was 
carried out with the High Energy Transmission Grating (HETG), which, in principle, gives 
a high-resolution source spectrum at higher energies. However, the dispersed spectra, 
integrated along the chip columns in the CC mode, are severely contaminated by the 
chip background, which strongly complicates their analysis. Therefore, we use only the 
zero-order image (on the back-illuminated ACIS chip S3) which provides a pulsar spectrum 
with resolution of about 0.1-0.2 keV. The archival pipe-line processed data are not corrected 
properly for the satellite wobbling (dither) and for the Scientific Instrument Module 
motion, which smears the 1-D image. We applied these corrections making use of the 



3 version of 31 Oct 2000; |http : //asc . harvard . edu/cal/Links/Letg/User/Hrc_QE/EA/ 
^ Preliminary results of this observation were presented by Stage et al. (2000). 
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algorithm kindly provided by Glenn Allen (see Zavlin et al. 2000 for details) and obtained 
one-dimensional (1-D) brightness distributions of the Vela pulsar and its PWN (bottom 
panel of Fig. 1). The upper panel of Figure 1 shows a 2-D image of this object obtained 
with the ACIS-S in the standard Timed Exposure mode (ObsID 128)0. This image is scaled 
and oriented in such a way that its horizontal axis is parallel to the 1-D image and its 
horizontal size is equal to the length of the 1-D distributions shown in the bottom panel. 
The 1-D profiles in different energy bands demonstrate that the pulsar is much brighter 
than the PWN background at < 1.5 keV, and it is discernible up to 8 keV, where the 
charged particle background takes over. 

For the analysis of the pulsar spectrum, we extracted counts from a 1-D 2'.'5-wide 
aperture. The nonuniform (in both the total brightness and spectrum) PWN background 
was found by interpolating between the adjacent regions on both sides of the pulsar. The 
estimated source count rate is 288 ± 4 ks~^, or 320 ± 5 ks"^ after correcting for the 90% 
fraction of point source counts in the 1-D aperture chosen. Because of the dither, the 
pulsar and the PWN were imaged on two nodes of the chip (node and node 1), which use 
different amplifiers and have different gains and responses. Therefore, we consider the two 
nodes as separate instruments while fitting the observed spectra. Since the ACIS response 
below 0.5 keV is still poorly known, we discarded the corresponding counts. 

We find that the ACIS spectrum of the pulsar does not fit with a one-component 
model, neither standard (BB, PL) nor a NS atmosphere spectrum, but it fits fairly well 
with two-component models (e.g., a PL plus a thermal component). 

The spectral parameters obtained from fitting the separate HRC-S/LETG and ACIS-S 

^Preliminary results of this observation were presented by Pavlov et al. (2000). See also 
tittp : / /chandra . harvard . edu/press/00_releases/press_060600vela . htnil| 
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spectra, as well as from the combined fit, are presented in Table 1. While fitting the LETG 
spectra, one cannot reliably find the parameters of the harder component, so we fix them at 
the values found from the ACIS-S fits. Similarly, the hydrogen column density riH cannot 
be determined reliably from the ACIS-S spectra with low-energy counts discarded, so we fix 
nn at the value found from the LETG fit. We see from Table 1 that the spectral parameters 
obtained from the separate fits are generally consistent with each other; the differences can 
be attributed to systematic errors caused by inaccuracies of the HRC/LETG and ACTS 
responses. The combined (LETG plus ACIS) fits (see Fig. 2) can be considered acceptable, 
in view of the current uncertainties of the instrument responses. Although we cannot 
completely exclude the possibility that the deviations from the continuous spectra might be 
due to some broad, shallow spectral features, the shortness of the exposures and the current 
knowledge of the instrument responses do not allow us to make definitive conclusions based 
on the data available. 

3. Discussion 

Thanks to the superb angular resolution of the CXO telescope, we are able to resolve 
the pulsar from its PWN and investigate its spectrum in the 0.25-8.0 keV range. Our 
analysis proves that the pulsar soft X-ray emission, at -E < 1.8 keV, is indeed dominated 
by a soft thermal component, so that the Three Musketeers (see Introduction) are now 
joined by D'Artagnan - the Vela pulsar. The parameters of the soft thermal component 
are drastically different for the BB fit^ {T{^ = 1.4-1.5 MK, = [1.9-2.4]d25o km — 
close to the values obtained by Ogelman et al. 1993 from the ROSAT spectrum) and the 

^ and are the temperature and the radius given by the standard BB fits; we use 
the superscript °° by analogy with the case of atmosphere fits, to emphasize that in both 
cases the quantities are given as measured by a distant observer. It should be noted that 
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hydrogen atmosphere fit {T^ = 0.65-0.71 MK, R°° = [14-17] ^250 km). For the BB model, 
one might assume that the observed radiation is emitted from small hot spots (e.g., polar 
caps heated by relativistic particles produced in the pulsar magnetoshere), whereas the R/d 
ratio obtained in the H atmosphere fit implies radiation emitted from the entire NS surface, 
with a lower temperature. 

Rigorously speaking, the surface (atmosphere) of a NS is not a black body, but a BB 
spectrum could mimic the spectrum of a heavy-element (e.g., iron) atmosphere if it is 
observed with low spectral resolution (Rajagopal, Romani, & Miller 1997). However, a lack 
of significant spectral lines in the LETG spectrum of the Vela pulsar hints that there are 
no heavy elements in the radiating atmosphere layers (although the exposure was too short 
to take full advantage of the high spectral resolution) . 

Because of the enormous NS gravity, the outer layers of the NS atmosphere should 
be comprised of the lightest element present. If, in the absence of H, the NS has a He 
atmosphere, we still expect spectral features in the observed range if the magnetic field 
is as strong as predicted by the standard radio pulsar models, S ~ 3 x 10^^ G. In such a 
field the He atmosphere is not completely ionized even at T ~ 1 MK, mainly because of 
the large increase of the ionization potentials (Ruder et al. 1994): / ~ 0.63 keV for the 
one-electron He ion (fa 0.5 keV with account for the gravitational redshift). The fact that 
we see neither a photoionization edge around E ~ 0.5 keV nor spectral lines at somewhat 
lower energies suggests that there is no He in the atmosphere. On the other hand, if even 
a small amount of H is present in the surface layers (~ 10^^ — 10^^ g of H over the NS 
surface would be enough to hide heavier elements completely), the hydrogen, with its lower 
ionization potential, / ~ 0.23 keV, should be strongly ionizied, and if even some neutral 

and R°° /d very weakly depend on the input parameters of the hydrogen atmosphere models 
(Zavlin, Pavlov, & Triimper 1998). 
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fraction is present, it will not show spectral features at £■ > 0.17 keV. Thus, the featureless 
spectrum we observe is consistent with the hypothesis of a hydrogen NS atmosphere. The 
observed part of the (continuum) spectrum of such an atmosphere decreases with E slower 
than the ideal Wien spectrum because the radiation at higher energies comes from hotter 
layers (the H opacity strongly decreases with frequency). As a result, BB fits to the H 
atmosphere spectra give a temperature higher than the true effective temperature and 
an area smaller than the true emitting area (e.g., Pavlov et al. 1995). If we adopt the H 
atmosphere hypothesis, the effective temperature of the Vela pulsar is below the predictions 
of the so-called standard models of the NS cooling (e.g., Tsuruta 1998), if even the pulsar 
is a factor of 2-3 older than its characteristic age (Lyne et al. 1996). 

An important result of our analysis is the detection of the hard spectral component. 
Because of the small number of high-energy photons detected, both the PL and thermal fits 
of this component are formally acceptable. However, the thermal fits yield an implausibly 
small size, ~ 10 m, and a very high temperature, ~ 10 MK, of the emitting region, and 
extrapolation of the model spectrum to lower and higher energies is inconsistent with 
optical and hard- X-ray observations. Therefore, we favor the nonthermal interpretation of 
the hard component — a PL (7 = 1.2-1.8 if the H atmosphere model is adopted for the soft 
component) which dominates at E >1.8 keV. With this interpretation, the overall X-ray 
spectrum of the the Vela pulsar is quite similar to those of the Three Musketeers, which 
show the nonthermal component dominating above 1.5-2.0 keV, with 7 ^ 1.3-1.6. It is 
interesting that the ratio of the nonthermal X-ray luminosity (estimated as if the radiation 
were isotropic, = Aird'^Fx) to the spin-down energy loss, L^/E = (1.7-2.6) x 10~^ in the 
0.2-8.0 keV range, or (2.1-3.4) x lO"'^ in the ROSAT range 0.1-2.4 keV, is much smaller 
than Lx/E ~ 10~^ for majority of X-ray detected pulsars (Becker & Triimper 1997). Thus, 
luckily for us, the magnetosphere of the Vela pulsar is very under luminous in the X-ray 
range — if its L^/E were as high as for other pulsars, the thermal radiation would be 
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undetectable in the phase-integrated spectrum. 

Since the Vela pulsar has been detected in the optical and gamma-ray ranges, where 
its radiation is certainly nonthermal, it is illuminating to compare its X-ray nonthermal 
spectrum with those in the other ranges. Figure 3 shows the energy flux spectrum of the 
pulsar from optical to gamma-rays. We see that extrapolation of the PL spectrum obtained 
in the atmosphere-|-PL fit matches fairly well with the optical and hard- X-ray (soft-7-ray) 
time-averaged fluxes. On the contrary, the slope of the PL component in the BB-I-PL flt 
is too steep to be consistent with the optical and gamma-ray data, which is an additional 
argument in favor of the H atmosphere model. Moreover, if we adopt the hypothesis that 
the spectrum has the same slope in the optical through X-ray range, we can include the 
optical points in the flt for the PL component, which constrains the photon index and the 
power-law normalization very tightly: 7 = 1.35 — 1.45, 

In the present paper we analyze only the phase-integrated spectra. Energy-integrated 
light curves, with at least three peaks per period, have been observed with the CXO HRC-I 
detector (Pavlov et al. 2000; Helfand, Gotthelf, & Halpern 2000). In fact, the phase- resolved 
spectral analysis, or energy-resolved timing, can provide much more information about the 
pulsar. For instance, we should expect that the pulsed fraction, pulse shapes (and, perhaps, 
even the number of peaks per period) are quite different at energies below and above 
1.8 keV because of the different emission mechanisms. However, because of the low pulsed 
fraction and complicated light curves, such an analysis would require an order of magnitude 
more pulsar photons than 16,000 photons collected in the described observations. 

We are grateful to Leisa Townsley and George Chartas for the useful advice on the 
analysis of ACIS data. Our special thanks are due to Glenn Allen, who provided the 
algorithm to correct the event times and coordinates in CC mode. We thank Gottfried 
Kanbach and Volker Schonfelder for providing the EGRET and COMPTEL fluxes in a 
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Fig. 1. — Bottom: One-dimensional profiles of the Vela pulsar and its PWN in different 
energy bands. One pixel is equal to 0'.'492. Top: 98" x 70" ACIS-S3 image of the Vela pulsar 
and the PWN, with the horizontal axis parallel to the orientation of the 1-D image. 

Fig. 2.— Combined fit of the HRC-S/LETG and ACIS-S3 count rate spectra with the 
two-component (magnetic hydrogen NS atmosphere plus power law) model. The dotted and 
dashed lines show the contributions of the thermal and nonthermal components. The best-fit 
parameters are given in Table 1. 

Fig. 3. — Multiwavelength (1 eV - 10 GeV) energy spectrum of the Vela pulsar. In addition 
to the CXO spectrum from this work, data shown include results from optical (Nasuti et 
al. 1997; U,B,V fluxes, without de-reddening), OSSE (Strickman et al. 1996), COMPTEL 
(Schonfelder et al. 2000) and EGRET (Kanbach et al. 1994) observations. The solid line 
shows a two-component (NS atmosphere plus power-law) fit to the observed CXO 0.25-8 
keV spectrum with spectral parameters — 0.68 MK, R°° — 15.5 ^250 km, 7 = 1.4, 
jV = 2.3 X 10~^, nH,2o = 3.2. The dotted fines correspond to same model spectrum corrected 
for the interstellar absorption and extrapolated to lower and higher energies. The dash-dot 
fines show the extrapolated optical and EUV spectra as absorbed by the interstellar matter. 
The thin solid lines demonstrate the uncertainty of the inferred power- law component. A 
typical error of the thermal component is shown by the vertical errorbar at the maximum of 
the spectrum. 
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Table 1: Parameters of the two-component spectral fits, with la uncertainties. 



Parameter 


HRC-S/LETG 


ACIS-S 


combined 




Blackbody + Power law 




T^OO 

-'bb 


1.4/ ± U.Ud 


1 1 _i_ n no 
l.Ol ± U.Uo 


1.49 ±0.04 


DOO 

-n-bb 


2.2 ± 0.3 


2.0 ± 0.2 


2.1 ±0.2 


7" OO 

-^bol 


i.D ± U.o 


1 /I 4- n Q 
i.4 ± u.o 


1.5 ±0.4 


7 


(2.2) 


1. 1 ± U.4 


2.7 ±0.4 


Af 


(6.0) 


n _i_ 1 Q 
D.U ± i.o 


10.7 ± 1.0 




(1.8) 


z.z ± u.o 


4.2 ±0.6 




i. / zh U.o 


(1.7) 


2.2 ±0.3 


xl [d.o.f.] 


i.u [loyj 


1 Q fl 1 Ql 

i.o [lloj 


1.1 [274] 




H atmosphere 


+ Power law 




T^OO 

-'eff 


U.D/ ± U.Uo 


U.D9 ± U.Uo 


0.68 ±0.03 


d 


2UU ± 25 


22U ± 2U 


210 ±20 


T OO 

-^bol 


3.9 ±0.3 


3.7 ±0.3 


3.8 ±0.3 


7 


(1.4) 


1.4 ±0.3 


1.5 ±0.3 


Af 


(2.3) 


2.3 ±0.7 


2.5 ±0.6 


Lx,pi 


(0.5) 


1.3 ±0.4 


1.5 ±0.4 




3.0 ±0.3 


(3.0) 


3.3 ±0.3 


xl [d.o.f.] 


0.7 [159] 


1.2 [113] 


1.0 [274] 



Note. — The values in parentheses were fixed during fitting. The temperatures 
and are in million kelvins (MK). The radius (in km) and the luminosities Lg^j and 
Lx,p\ (in 10^^ erg s~^) are related to d — 250 pc. The nonthermal luminosities L^^pi are 
in the bands 0.2-2, 0.5-8, and 0.2-8 keV for the LETG, ACIS-S, and combined spectra, 
respectively. 7 and A'' are the photon index and normalization constant of the power-law 
component: dN/dE — J\fE~'^, E in keV, N in 10~^ photons/(cm^ s keV). The hydrogen 
atmosphere models were calculated for S = 3 x 10^^ G, M = I.4M0, = 10 km {g^ = 0.766, 
= 13.05 km). 
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